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^* Abstract 

*7' A complete error analysis of variational integrators is obtained, by 

blowing up the discrete variational principles, all of which have a singu- 

i-Q larity at zero time-step. Divisions by the time step lead to an order that 

is one less than observed in simulations, a deficit that is repaired with the 
ri help of a new past-future symmetry. 

co 1 Introduction 

> 

^^•O We consider a regular Lagrangian system L: TQ — > M, and associated Euler- 

Lagrange equations 

P"- dtdq 1 dq l ~ 

o 

00 Constraints if present arc holonomic and are incorporated into the configuration 

(^ manifold Q. Standard numerical integrators are insensitive to the specialties of 

£>. such conservative systems because they discretize Euler-Lagrange equations as 

? *~j they do any other differential equations. 
rN Variational integrators [7J|8j[TT] rather discretize Hamilton's variational prin- 

03 ci P le 

b 

L\q'[t)) dt = 0, q(a) and q(b) constant 

(q'(t) £ TQ includes both coordinates q l and v l in this notation). Such dis- 
cretizations can be obtained by fixing a time step h > and approximating 
the integral with a finite sum, resulting in a finite dimensional constrained op- 
timization problem — a discrete Hamilton's principle. The function under the 
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sum, called the discrete Lagrangian, is an approximation to the type 1 generat- 
ing function 

S h (q + ,q-)= [ l L{F t (A h (q+,q-)))dt 
Jo 

where Ah(q + ,q) is the initial velocity at q~ that arrives at q + rs q~ after 
time h, and F t is the flow of the Euler-Lagrange equations. The critical points of 
the discrete Hamilton's principle satisfy discrete versions of the Euler-Lagrange 
equations, which then define a numerical integrator for the original Lagrangian 
system. 

Marsden and West [8] consider the local existence and uniqueness and the 
local error analysis (order) of variational integrators. But these problems are 
quite subtle: the map A^(g + ,g _ ) is singular as h — ► + , since arbitrarily large 
velocities are required to go from fixed q~ to fixed q + in vanishingly small time. 
We show (Section |2| that terms not controlled by definitions of the order of 
a discrete Lagrangian system, enter in the discrete Euler-Lagrange equations, 
such that a complete error analysis requires additional argumentation beyond 
that in [3]. 

We obtain a complete local existence and uniqueness, and local error anal- 
ysis, for variational integrators of Lagrangian systems L: TQ — ► R. To system- 
atically treat the singularity, we use the discretizations of Cuell and Patrick [5] , 
which 

1. discretize the tangent bundle of Q using finite curve segments; and 

2. place the discrete variational principle in the velocity phase space TQ. 

This shifts the singularity, to a degeneracy in the map that sends the curve 
segments to their endpoints. Viewing discrete tangent vectors as curve segments 
is generally consistent with viewing discretizations in general as attaching to a 
manifold finite rather than infinitesimal objects [TJ [2J. Our discrete tangent 
bundles are similar to the groupoid based constructions of discrete phase spaces 
for Lagrangian systems in |10j . 



To prove existence and uniqueness (Theorem 3.7 1, we blow up the varia- 
tional principle at h = 0, converting it to a smooth perturbation of the trivial, 
nonsingular optimization problem, with objective L(v q ) + L(w q ) and constraint 
v q + w q = constant. Divisions by h lead the arguments to an order that is one 
less than observed in simulations. But the past and the future occur symmetri- 
cally in the blown-up variational principle, giving a new Z2 symmetry (exchange 
of v q and w q ), from which we prove the observed order follows by a nontrivial 



cancellation (Theorem 4.7). A similar cancellation occurs in the error analy- 
sis in [5], but that development is restricted to type 1 generating functions of 
kinctic-plus-potential systems, and it is framed non-variationally, on cotangent 
bundles. 

The use of extended phase spaces, particularly the Hamilton-Pontryagin 
variational principles, where there is also additional kinematic freedom, is an 



emergent interest [31 [T2J Q2]. In this work, the extra kinematic freedom in 
the discrete tangent bundle phase space, as opposed to configuration space, is 
necessary as a place for the discrete Lagrangian system to go, in the limit as 

We obtain semi-global results: Theorem |3.7| asserts a well defined discrete 
evolution for arbitrarily high velocities, if h is sufficiently small i.e. the discrete 
evolution is defined for (h,v) in an open neighborhood of {0} x TQ. We an- 
chor our arguments directly to the discrete variational principles, avoiding use 
of auxiliary constructions, such as the discrete Euler-Lagrange equations, or 
canonical contexts on the cotangent bundles. These are important and useful, 
but it is also important and useful to obtain the discrete variational picture 
stand-alone, as a coherent whole. For all these purposes, and also to make clear 
the geometry underlying the error analysis, our theory is coordinate free. 

We extensively use the development in [4 . In this paper, all manifolds are 
assumed paracompact and Hausdorff. 

2 Preview by example 

The following example-oriented preview provides a good understanding of our 
motivations, and of the results of the general theory in Sections |3j and [4] 

2.1 Definition 

We consider Q = E — {q}, with TQ = M 2 = { (q, v) } and Lagrangian 

L = l -g(q)v 2 - V(q) (2.1) 

where g(q) > and V(q) are smooth functions. Such systems describe a particle 
moving in Euclidean space under the influence of a potential, and constrained 
to some curve parameterized by the variable q. 

The discrete phase space for this system is Q x Q = { (q + , q~ ) } , and a 
discrete Lagrangian may be defined by 

L h (q+,q-) = hL (q- t 9+ ~ q ) + ah(q+ - q'), (2.2) 

where a is a constant that is necessary for the purpose of the example. The 
discrete action on a sequence qo , gi , q-i € Q is 

Sh = L h (q 1 ,q ) + L h (q 2 ,q 1 ). (2.3) 

From this data, a variational integrator is defined: given an initial discrete 
state (qi,qo), a time h advanced discrete state (q2,qi) is computed from the 
requirement that (qo,qi,q2) is a critical point of Sh, under the constraints that 
qo and q 2 are not varied. This is the standard setup of discrete mechanics [51ITTJ. 



2.2 Order 

Substituting into Lh the flow of Euler-Lagrange equations i.e. 

q~ = q, q + =q + vt + 0(t 2 ) 

gives, after putting t = h, 

L h (q, v) = hL(q, v) + ah 2 v = L h (q, v) = hL(q, v) + 0{h 2 ). (2.4) 

The integral of L along the flow is 
r h 

/ L(q + 0(t),v + 0(t))dt=hL(q,v) + 0(h 2 ). (2.5) 

Jo 

By standard definition, a numerical integrator yi+i = F(yi,t) of differential 
equation y' — f(y, t) is order r if the local truncation error y± — y(t + h) is 
0(h r+1 ) as h — > 0. According to Equation 2.3.1 of [8 , L/j is by definition order 1 



since ( |2.4| and (2.5) agree through order ft, 1 . So, Item (1) of Theorem 2.3.1 of [8 
asserts that the variational integrator defined by Lh should be order 1. 

Indeed, the critical points of Sh are exactly the solutions of the discrete 
Euler-Lagrange equations 

F-L h (q 2 ,q 1 ) = F + L h (q 1 ,q ) (2.6) 

where 

jp+r - dLh jp-r - dLh 

Oq + Oq~ 

are the discrete Legendre transforms. One computes 

F+L h = pj{q~){q + -q~) + axh, 

F~L h = ^s(<T)(g + - q~) + hV'(q-) +axh- l- g \ q -)( q + - q ~f, 



after which Equations (2.6) become 



-r9(qi)(Q2 - qi) + hV'(qi) +axh- ^g'(qi)(q 2 - qif 



T9(qo)(qi -qo)+axh. 



(2.7) 



Assuming qi — qo = 0(h), this implies 

which is evidently order 1, after comparison with the Euler-Lagrange equations 



of L defined by (2.1 ), which are 



" 1 " l f V'(q) + y(qy 



dP g(q) V W 2* 

Item (1) in Theorem 2.3.1 of |8] asserts that the discrete evolution should be 
order 1, and it is. 



2.3 Discrete Legendre transforms 

Something, however, is not quite right: the terms axh in the discrete Lagrangian 
Lh are not controlled by the definition of an order 1 discrete Lagrangian, be- 



cause they enter Equations ( 2.4 1 and (2.51 only at order 2. In the discrete 



Euler Lagrange equations (2.7 'I such terms appeared twice, and at the same or- 
der as the potential, which they could have significantly altered. They appeared 
symmetrically on both sides of the equation, and therefore did not affect the con- 
sistency of the variational integrator. Nevertheless, they did affect the separate 
terms appearing as discrete Legendre transforms in the discrete Euler-Lagrange 
equations. 

Item 2 in Theorem 2.3.1 of [S] asserts that the discrete Legendre transforms 
F + Lh and F~Lh should, after the substitutions q = q~ , q = q~ + hv, be or- 
der 1 consistent with the exact Legendre transform of L. This shifted the term 
ah(q + — q~) to order 2, because it enforced q + — q~ = 0(h). So that term 
did not affect the order 1 consistency of the discrete Lagrangians. However, 
that term did affect the discrete Legendre transforms at order 1, because the 
derivatives in the Legendre transform removed q + —q~ . Therefore the definition 
of an order 1 discrete Lagrangian did not control the order 1 consistency of the 
discrete Legendre transforms. The fault may be assigned to the \/h singularity: 
It was decided to clear that singularity in the definition of the order of discrete 
Lagrangians, because the definition of order is necessarily about h = 0. When 
the singularity was cleared, the term ah(q + — q~) was knocked out of the defi- 
nition, but it affected the discrete equations of motion. The situation was saved 
only because ah(q + — q~) appeared in the discrete equations of motion as a pair 
of identical twins. 

Suppose we view variational integrators just as the discrete Euler-Lagrange 
equations, discarding the variation theory altogether. Possibly for some pur- 
pose of numerical implementation, imagine being motivated to use different 
discrete Lagrangians in the separate terms of the discrete Euler-Lagrange equa- 



tions i.e. to use rather than (2.6), the equations 



F L h {q 2 ,qi) = F + L h { qi ,q Q ), 

where L ^ L. Then the order of the corresponding numerical integrator would 
almost certainly be reduced by 1, despite the fact that each Lagrangian would 
separately satisfy all prescribed order conditions. And, in the case of order 1 
Lagrangians, the result would almost certainly be an inconsistent integrator. 

2.4 Example analysis: blow up 

To analyze the discrete variational principle as h — > 0, change to new variables 
q,v,q, v defined by 

q = q , v= — : , q = qi, v= — : . (2.8) 



There are four new variables but originally there were only qo,qi,q 2 , so there 
will be a new constraint. Defining 

9h(l>v)=Q> d£(q,v) = q + hv, 



one has, from (2.8 1, 

9u(Q,v) = Qo, d£{q,v) = qi=d^(q,v), d h h (q,v) = q 2 . (2.9) 

The new constraint is the middle of these; the variations in the discrete varia- 
tional principle are restricted to annihilate the derivatives of the outer two. The 
action becomes 



Sh = L h (q,v) + L h (q,v), L h (q,v) = hL(q,v) + ah 2 v, 



(2.10) 



which, together with the constraints (2.9), is an equivalent variational principle 
onTQxTQ = {(q,v),(q,v)}. 



We will desingularized the variational principle defined by (2.9 1, (2.101, in 
three stages: 

1. The critical points are unchanged by multiplication of Sh by any constant. 



So Lh in (2.101 can be replaced with its division by h. This results in a 



replacement of Sh with its division by h, although the formula for Sh is 
unchanged. The variational principle, objective and constraints, becomes 



Sh = L h (q, v) + L h (q, v), L h (q, v) = L(q, v) + ahv, 
9h(q,v)=q , d£{q,v) = q r = d^(q,v), d£(q,v)=q 2 . 



(2.11) 



2. The middle constraint of (2.111 is q + hv = q. This constraint can be 



imposed by substituting q, resulting in a principle on the variables (q, v, v). 
The variational principle becomes 



Sh = L h {q, v) + L h (q + hv, v), L h (q, v) = L(q, v) + ahv, 
q = qo, q + hv+ hv = q 2 . 



(2.12) 



3. At h = 0, the fixed endpoint constraints of (2.121 degenerate to the same 



function, namely (q,v,v) i— ► q. To remove this degeneracy, these can be 
post-composed by any smooth bijection, such as 



9o + 92 



92 -Qo 



The variational principle becomes 

Sh = L h (q, v) + L h (q + hv, v), L h (q, v) = L(q, v) + ahv, 

h, „. 

q+ -(v + v) = q, v + v — z. 



(2.13) 



The desingularization is complete because the principle (2.131 is smooth and, 
as will be seen, nonsingular, through h = 0. The blown-up variational principle 



is obtained from (2.131 by substituting h — i.e. 



Sq — L(q,v) + L(q,v), v + v = z. 



(2.14) 



perturbation of (2.141 



For small h, the variational principle (2.131 may be regarded as a continuous 



2.5 Example analysis: discrete existence and uniqueness 



Re-introducing the details of the example Lagrangian (2.1 1, the blown-up vari- 
ational principle becomes 



S = 9 ^(v 2 + V 2 )-V(q), 



V + V = z. 



(2.15) 



Substituting v — z — v and differentiating gives, since q is constant in (2.15) 



dS 
dv 



= g(q)(2v-z) = 0, 



so, given z, there is the solution v 
critical point is 



z/2. The second derivative at this 



d 2 S 

dv 2 



2g(q), 



which is positive, so the critical point is nondegenerate. Nondegenerate critical 
points vary smoothly with parameters and constraint values. In particular, for 
any go an d zq, the original principle (2.9 1, 
(v,v) 



(2.101 has a unique critical point 
go and z 



7(/i, q, z) for all sufficiently small h, and all q « go an d z ~ z o- 7 j s 
smooth and at h = its graph is an open set of the diagonal of TQ x TQ i.e. the 
graph of the identity map. So, for sufficiently small h the graph of 7, and hence 
the variational principle, locally defines a near- identity time-advance map of 



TQ to itself i.e. a discrete evolution. For the general results, see Definitions 3.1 



and 3.3 and Theorem 3.7 



2.6 Example analysis: order 

Error analyses of variational integrators has purpose to give conditions on the 
basic data of the discrete variational principles that are sufficient to imply that 
the corresponding numerical integrators are order r. This follows if: 

1. exact discrete variational principles are obtained that give the exact solu- 
tion of the continuous Lagrangian system; and 



2. two variational principles that are order r 
that are order r + 1 consistent. 



1 consistent have solutions 



A discrete variational principle is then defined order r + 1 if its objective and 
constraints are order r + 1 consistent with those of the exact discrete variational 
principle, and then Item (2) implies that the corresponding numerical integrator 
is order r. 



Item (1) seems not to be clarified by specializing to Example (2.2), but 
rather seems more clear in the geometric setting; see Lemma |4.2| and Theo- 
rem [43] Essentially, in gross aspect, solutions of an exact discrete variational 
principle correspond to cornered solutions of the continuous variational prin- 
ciple. However, by the well known Weierstrass-Erdman conditions [6], regular 
variational principles do not have corners i.e. critical points are unaffected if 
the space of solutions is enlarged to allow corners. Therefore, solutions of the 
continuous variational principle are solutions of the exact discrete variational 
principle. But then the critical points of both principles are the same because 
of uniqueness i.e. the solutions of the discrete principle are exact. 

Item (2) is clarified by the example context. The desingularized princi- 



ple (2.13) has discrete action 

Sh = L(q, v) + L{q + hv, v) + a x hv + a x hv. 

At a solution h = and v = v, Sh is symmetric under the exchange v <-> v, not 
just at order hP, but also at order h 1 , because the original minus the exchanged 

is 

S(q,v,v) -S(q,v,v) 

= L{q, v) + L(q + hv, v) — L(q, v) — L{q + hv, v) + ah{v — v) 

dL dL 

= L(q,v) — L(q,v) + h — (q,v)v - h — (q,v)v + ah(v — v) + 0(h 2 ), 



which is 0(h 2 ) at v — v. So the constraint q = q+hv of (2.9), and the term axh, 
affect the solutions of the discrete variational principles only symmetrically at 
order h . Considering the solution of the discrete principle as a perturbation 
from h = i.e. a perturbation of the graph of the identity map, the effect on 
the solutions is symmetric to order h 1 . But such a symmetric alteration does 
not change, at the same order, the function defined by the graph. In particular, 
the solutions of the variational principle are unaffected, at order h , by the term 
axh. For the general results, see Definitions |4.1| and |4.6[ and Theorem |4.8| 

3 Discrete existence and uniqueness 

Our discretizations of Lagrangian systems depend on discretizations of tangent 
bundles of manifolds Ai, by assignment of curve segments in M. to tangent 
vectors of M. [S] . We require a parameter h such that TA4 is obtained in the 
limit h — > + , so we posit a map ip(h, t, v m ), with values in Ai, and obtain the 
curve segments t i— > ?p(h,t,v m ). The definition builds in some flexibility. The 
curve segments are generated as the variable t ranges over intervals of length h, 



with otherwise unrestrained endpoints: for example, [0, h] and [— h/2, h/2] are 
common choices which arc both accommodated. 



Definition 3.1. A C k discretization of TM, k > 1, is a tuple (4>,a + ,a ), 
xTM^M, a + : [0,o) -^ R> , a~:[0,a) 



where 

V>: 17 C 

are such that 

1. ip is continuous, U is open, and {0} x {0} x TM C 17; 

2. a + ,a~ are C , and a + {h) — a~(h) = h; 

■-) ,; 

3. ip(h, 0, v m ) = m, and ^-(/i, 0, w m ) = u m ; 

ot 



d'l (v m ) = ip(h,a + (h),v n 



>-<o, 



. the boundary 


maps defined by 




dh( v m) = ip(h,a~(h),v m ), 


d+(v 


are C k in (h, v m ) and 




d 
dh 


d£(v m ) = a + v m , 

h=0 


d 
dh 


where 






■ + da+ ln\ 
a = ^ (0) ' 


6r 



d h (v„ 



a v„ 



(3.1) 



(3.2) 



h=0 



da 
dh 



(0). 



Remark 3.2. Putting h = in a + (/i) — q (h) = h gives a + (0) — a (0) = 
because a + > and a~ < 0. If ip is a C 1 map in all its variables then Equa- 



0,0+ 



tions (|3.2p are superfluous because they follow by differentiating Equations (3.1 1. 

~ = tq, and differentiating a + {h) — a~{h) = h 

i. n 



Also, note that at h 
at h — gives a + — a 

Given a discretization of the tangent bundle of configuration space, one only 
need add an appropriate discrete Lagrangian to obtain a discretization of a La- 
grangian system. See [5] for more explanation, and a general development along 
the lines used in this paper, of discrete Lagrangian mechanics and discretiza- 
tions of Lagrangian systems, extending to nonholonomic systems. Here we do 
not require the full generality, and what is required, specialized to the holonomic 
case, is collected in Definition 13.31 



Definition 3.3. A C k discretization, k > 1, of a Lagrangian system L: TQ 
M., is a tuple (L^, tp, a + , a~) where 

1. (i/j,a + ,a~) is a C k discretization ofTQ; and 



2. L h : TQ 



is C k in (h,v q ) is such that Lh(v q ) — hL{v q ) + 0(h 2 ). 



(h, Vq, vo) elx TQ x TQ is critical if v = Vq, v = vo is a critical point of the 
variational principle 

j S h = L h (v) +Lh(v), 

[ dr(v) and d£(v) constant, and d£(v) = d7{v). 

A discrete evolution is a map F defined on an open subset of R X TQ such that 
(h, v, F(h, v)) is critical for all (h, v) in the domain of F. 

Remark 3.4. Precisely, (h,vo,vo) is critical if 

dSh(v,v)(Sv,5v) = and d^(v) — d7(v) 

for all Sv eT v Q and 5v £ T% Q such that 

Td^(Sv) = 0, Td+(Si) = 0, Td+(Sv) = Td^(Sv). 

We have placed the discrete variational principle in sequences in velocity phase 
space, so there is a discrete analogue of the (continuous) first order constraint 
q(t)' = v{t) where q(t) = igo v(t) i.e. the successive curve segments must 
join to make a continuous whole. The fixed endpoint constraints correspond 
to d~^(v) and d£(v) constant, which affect only the variations; the actual con- 
straint values are not specified. This provides, as is usual in both discrete and 
continuous Lagrangian mechanics, the necessary freedom to accommodate ini- 
tial conditions. □ 

Remark 3.5. Suppose that L® xC (fi + , q~) is a discrete Lagrangian as defined 



in [8] and let ip be a discretization of Q as in Definition 3.3 As shown in [5], 
^fh(v) = (d£ (v) , d^ (v)) is a diffeomorphism from a neighborhood of the zero 
section of TQ to a neighborhood of the diagonal of Q x Q. Lh = L^ ° ^h 
is a discretization as in Definition |3.3| and conversely, any such discretization 
defines a discrete Lagrangian as in [5] by L^ xQ = L h o ^^ l . In any case, the 
discrete mechanics on TQ and QxQ are equivalent because they are conjugated 
by iff. For more details, see Section [5| □ 

At h — 0, the discrete action is singular, and the constraints are degenerate, 
because 

1. at h = 0, L h (v) = 0, so Sh(v, v) = L h (v) + L h (v) — 0; and 

2. at h = 0, dt{v) = d7{v) is tq(v) — tq(v), and on that constraint the 
fixed endpoint constraints d7 {v) = q~ and dt(v) = q + are replicates. 

The necessary blow-ups rely on a technical result of [1] that we recall as Propo- 



sition 3.6 below. This is an invariant version of the elementary calculus fact, a 
version of L'Hospital's rule, that if f(x) = f(x)/h(x) where f(x) and h(x) are 
C 1 functions such that /(0) = and h'(0) ^ 0, then f(x) can be continuously 
extended through x = by defining /(0) = f'(0)/h'(0). 



10 



If 7r : E — > Ai is a vector bundle, and z S T$ m E, then we denote the horizon- 
tal and vertical parts of z by horz £ T m M and vertz e E m , respectively. We 



denote the zero section of E by 0(E). Also, the statement of Proposition 3.6 
uses the convention that a pair (Ai,hj^) is called a manifold when .M is a 
manifold and /i_vi : M. — > R is a submersion. 

Proposition 3.6. Let (A4, /ix) ai3< ^ M be a manifolds, and let it: E — > TV 
be a vector bundle. Suppose that f : U £ Ai — > E is C k , k > 1, and that 
/(m) € O(-E) whenever /ijk(to) = 0. Then for all m such that /^(m) = 0, 
there is a unique e(m) £ E n /f^ m ^ such that 

vertT„J(v m ) = (dh M (m)v m )e(m), v m £ T m M. 

Moreover, the function f : U —* E defined by 

i r — 7— t, ^M(m) ^0, 
/(m) ee { /ix(m) 

e(m), /ix( TO ) = 0? 

isC*" 1 . 

Theorem |3.7| is the main result on discrete existence and uniqueness. This 
goes beyond just finding regularity conditions on the discrete Lagrangian that 
imply existence and uniqueness, because it only relies on regularity of the contin- 
uous Lagrangian — it is an analysis of the limit h — > 0. For the discrete analogue 
of regularity in our context, and the corresponding existence and uniqueness 
result, see [5], and in the standard context that is Theorem 1.5.1 of [8j. 

Recall that the first and second fiber derivatives of L : TQ — ► R are FL(v q ) = 
D(L\T q Q)(v q ) and F 2 L(v q ) = D 2 (L\T q Q)(v q ), and that L is called regular if 
F 2 L has values only in the nondegenerate quadratic forms on the fibers of TQ. 
If A is a set, then we denote by A(A x A) the diagonal of A x A. 

Theorem 3.7. Let (Lh,i/j,a + ,a~) be a C k discretization of a regular La- 
grangian system L: TQ — > R, k > 2. Then there are neighborhoods W C 
RxTQof{0} xTQandU £RxTQxTQof{0} xA(TQxTQ) such that, 
for all (h,v) £ W, h > 0, there is a unique v £ TQ such that (h,v,v) £ U and 
(h,v,v) is critical. Moreover, U and W can be chosen such that F: W — » TQ 



dchned by 



F(h,v) 




isC 



fc-i 



Proof. The blow-up of Lh is immediate: set 

1 



L(h,v q ) = { h LhM > h *°> 
L(v g ), h = 0. 

11 



L is C k 1 by Proposition 3.6 



The constraints are blown-up by imposing d^(v) = d^ (v), after which d^(v) 
and dfr(v) are O(h) close and their difference can be usefully divided by h. For 
this, observe that both d£ and d^ are submersions on TQ when h = 0, since 
do( v q) = 9 an d do( v q) = 9' so there is a neighborhood i3 {0} x TQ on which 
both d£ and 9^~ are submersions. Consequently 

C = {(h,v,v):(h,v)€A,(h,v)€A,d+(v) = d^(v)} 

is a submanifold of M x TQ x TQ. Also, there is a tubular neighborhood 

C: W ^ C £ -> W QxQ C Q x Q = { (g+, <T) } 

of the normal bundle E = {(v q , —v q ) : v q 6 TQ] to the diagonal A(Q x Q) of 
QxQ, which satisfies 

vertrr^+.i;-) = Q(»+ - «-), i(v" - <)J . (3.3) 

The purpose of C _1 is to compute the difference between two nearby elements 
of Q. For example, if Q = M" we can use ((v q ,—v q ) = (q, (?) + (w g ,— w g ), and 
then 



i/„ + _^ ... .^ .. / 9' +? ? 



i.e. the fiber part of E corresponds to the difference. Just below, in the definition 
of <p, scalar multiplication of E by 1/h will be used to blow up the difference. See 



the proof Proposition 1.9 of [5] for more details about arranging Equation (3.3 1. 
Define ip:C^Rx£by 

<p(h,v,v) = { , 

h, -(v + v, —v — v) ] , h = 0. 

and define yj and /i£ on C by 

ip(h,v,v) = C _1 (5jJ"(ti),5^ («)), h c {h,v,v) = h. 

If (h,v,v) S C and /i = then tq(u) = 9q ( v ) = ^o"(") = T s(^)- Hence for all 



(h,v,v) £ C, ip{h,v, v) = if he{h,v,v) = 0. By Proposition 3.6 for all u,i) 
there is a unique e(w, v) such that 

vert T<^(0, u, v)(Sh, 8v,Sv) = [e(v ) v) ) —e(v ) v))dhc(0 > v,v)(Sh,Sv,dv) 

= (e(v,v),-e(v,v))5h, 
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for all (5h,Sv,Sv) G T^ 0v ^C. By Items [3] and PA of Definition 3.1 (Sh,5v,Sv) G 
T(o,v.v)C if and only if 

Ttq(Sv) + dha+v = Ttq(5H) + 5ha~v (3.4) 



and using Equation (3.3), and the definition of ip, 



vertTp(0,v,v)(Sh,Sv,Sv) — -(w,—w), 
where 

w = Ttq(5v) + 5hd + v — Ttq(Sv) — bhcTv. 
It follows that e(y,v) can be found by imposing 

-(Ttq(8v) + Sha + v - Ttq(Sv) - ShoTvj = She(v,v) (3.5) 



for all 6h, 5v, and Sv which satisfy Equation (3.4). Using A3 At to replace 



Ttq(Sv) — Ttq(Sv) in Equation ( |3.5| gives 
2 1 



- 1 5ha + v — 5ha v + 5ha + v — 5ha v) — -6h (v + v) = 5he(v,v) 



3.6 



so e(v, v) = 2 (v+v), after which the definition of ip at h — 0, and Proposition 
applied to tp, imply ip is C fc_1 . 

If h > 0, then finding the critical points of the discrete variational principle 
is equivalent to finding the critical points of L\(p~ 1 (h,z q ,—z q ). If h = 0, then 
the latter is the problem of finding the critical points (vq,vo) of L(v) + L(v), 
v, v G T q Q subject to the constraint |(t> + v) = z q . These are the v and 
v — 2z q — v such that L(v) + L(2z q — v) has a critical point at v i.e. such that 

FL(v) - FL(v) =0, v = 2z q - v. 

This has the solution v = v = z q , at which the Hessian of L(v) + L(2z q — v) 
is 2F 2 L(v), which is nondegenerate. Thus there is a manifold of nondegen- 
erate critical points parametrized by z q G TQ. Semiglobal persistence of 
these critical points follows by Theorem 2 of jl] i.e. there are neighborhoods 
U D {0} x A(TQ x TQ) and V 2 {0} x E, and a C*" 1 map 7: T> -> L>, such 
that for all (h, z qi —z q ) G V, j(h, z q ,—z q ) is the unique critical point in [/ of 
£|0 -1 (/i, «„-«,). 

At ft- = 0, ^(h, z q , ~z q ) = (0,z q ,Zq), and the image of 7 forms the graph 
of the identity map of TQ. Consequently, for small h, 7 determines a map F 
because 7 has image a graph. The technical statements in the Theorem to this 
effect are immediate from Proposition 5 of 0], applied to the map 7r 2 3 07, where 
V2a(h,v,v) = (v,v). □ 



Remark 3.8. The proof of Theorem |3.7| shows the blow-up at ft. = of the 
discrete variational principles gives the variational principles with action L(v) + 
L(v), where v and v are constrained (1) to be in the same fiber of TQ, and 
(2) such that v + v is constant. The blown-up variational principle is past-future 
symmetric i.e. symmetric under the exchange of v and v. □ 
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4 Order 

If the curve segments of the discretizations of the tangent bundle of Q are ob- 
tained from the base integral curves of the Euler-Lagrange vector field Xe, 
and the discrete Lagrangian is the classical action, then we obtain the exact 
discretizations of Marsden and West |5]. Exact discretizations are important 
because they exactly generate the flow of the Euler-Lagrange equation of the 
continuous Lagrangian, (Theorem 1.6.4 of [5]), so that the order of a discretiza- 
tion can be controlled by reference to its order of consistency with an exact 
discretization. 

Definition 4.1. An exact discretization of a Lagrangian system L: TQ — > R 
is a discretization (L/j, tp, a + , or) where ip and Lh satisfy 

1. ip(h,t,v q ) = T Q (F t XB (v q )), where Ff E is the Bow ofX E ; and 



~dt 



2. L h (v q ) = / Lo —(h,t,v q )dt. 



(h) 

Understanding of exact discretizations can be anchored to the variational 
principles. For an exact discretization, curve segment discretizations lead to a 
picture of piecewise smooth solutions built of segments of the exact flow. The 
difference between an exact discrete and continuous variational principles is that 
the discrete principle allows for piecewise smooth solutions — it is, after all, con- 
structed from segments. But the well known Weierstrass-Erdman conditions [6. 
imply that additional critical curves are not obtained by allowing corners. For 
example, corners do not occur in Riemannian geodesies because they are trian- 
gles for which there would be a locally shorter path along the hypotenuse than 
two of the sides. Lemma |4.2| and Theorem |4.3| provide the formal statements 
and proofs along these variational lines cf. Theorem 1.6.4 of |8]. 

Lemma 4.2. Let (Lh,ip,a + ,a~) be a C k exact discretization of the C k La- 
grangian system L: TQ — > R, k > 2, suppose that the integral curve of 
Xe through v € TQ is dchncd for times in [a~(h),a~(h) + 2h), and set 
v = F^+( h) (v). Then (h, v, v) is critical. 

Proof. The variational derivative of the action 

S= I L(q'(t))dt 

J a 

can be written [5j [7] as 

dS(q(t))-5q(t) = J 6L(^Ysqdt + 



Fl ( *],>«,/, 



(4.1) 



where SL is locally 



i d dL dL ' 

' )L= ^-dtW + W ]d9 
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Assuming Sv and Sv' satisfy the constraints 

Tdk(Sv) = 0, Td£(6v)=0, Td+{5v) = Tdl{8v), (4.2) 



and remembering that SL = along a solution, one applies Equation (4.1l to 
each of the integrals in 

pa+(h) j-a + (h) 

S h (v,v)= Lo(F t XE (v))'dt+ Lo(F t XE (i))'dt, 

Ja-(h) Ja-(h) 



obtaining, from Equation (4.2) and F^+r h \(v) — F Xb , }i (v), 



dS h (v,v)(Sv,Sv) = FL{F*f {h) (v)) Td+{Sv) - FL(F** (h) (v)) Td^(6v) 

= 0. □ 

Theorem 4.3. Let (Lh,ip,a + ,a~) be an exact discretization of a regular La- 
grangian system L : TQ — > E. Then there is a neighborhood U C E x TQ x TQ 
of {0} x A(TQ x TQ) such that, for all (h, v,v)&U with h>0,i = F Xe (v) 
if and only if (h, v, v) is critical. 

Proof. Let W C E x TQ and U C E x TQ x TQ be as in the statement of 
Theorem |3.7| Possibly by shrinking W, one can assume that 



1. for all (h,v) G W, the integral curve of Xe through v is defined for times 
in [a~(h), a~(h) + 2h]; and 



2. (h, v, FQ E {v)) G U for all (h, v) G W. 



4.2 



Con- 



If (h,v,v) G f and v = F h E (v) then (h,v,v) is critical by Lemma 

versely, if (h,v,v) is critical then so is (h,v,F h E (i>)), so w = F^ E (v) by the 

uniqueness in Theorem |3.7| □ 



We will avoid coordinates in order computations on manifolds by using the 
calculus of residuals developed in |3]. Suppose (M,fiM) an d A/" are manifolds. 
If fii M. — » A/", i = 1,2, are such that /i = /a on /i^"(0), then define /2 = 
/i + 0(/i^ / i), r > 1 if, for all mo G ft^(0), there is achart v at no = fi(mo) G A", 
and there is a function (^/) y defined near m , and continuous at m , such that 

v {h( m )) - "(AN) = h M {m) r (5f) v {m), 

for all to in some neighborhood of mo. As is easily shown, (Sf)v(mo) trans- 
forms as a tangent vector as v is varied, and therefore defines an element 
res r (/2, /i)(mo) G Tft mo \J\f, called the residual. Proposition 4.4 which is a 
specialization of Proposition 3 of [3] is the key result used to compute residuals 
without the invocation of local charts. 
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Proposition 4.4. Let (M,1im), (A/ - , hj^f), and V be manifolds, and suppose 
fi\ M. — > N and gi\ M — > V. i = 1,2 are C 1 and satisfy hj^ o / 4 = /i^(, 
/ 2 = /i + O(fcJU), and 52 = 3l + 0(/fr). Tlien 52 o / 2 = ffl o f x + O(feJU). 
Moreover, if hj^{m) = and n = fi(rn), then 

res r (g 2 o/2,ffi°/i)M = res r (g 2 ,gi)(n)+T n5 i res r (/ 2 ,/i)(m). 

As has been stated, the central issue is a decrease in the order of accuracy, 



essentially due to a division by h. Proposition 4.4 which is yet another result 



of [4], tracks this in the context of Proposition 3.6 



Proposition 4.5. Let (Ai, /i») and N be a manifold s, le t tv: E — > J\f a vector 
bundle, and suppose fi and fi are as in Proposition 3.6 with k > r. Then 
h = A + 0{h r M X ) if h = h + 0(h r M ), r > 2. Moreover, res r (/ 2 , A) takes 
values in the vertical bundle of E and res r_1 (/ 2 , /i) = res r (/ 2 , /i). 



We come to the main objective, accuracy, which is the order to which an 
evolution map of a given discretization of a Lagrangian system agrees with 
the continuous flow of that Lagrangian system. This can be approached by 
analyzing the order that two discretizations agree, since, by Theorem |4.3[ the 
continuous flow is obtained from the exact discretizations. Together with the 
equivalences of Section [5j this suffices to repair the proof of (3) implies (1) in 
Theorem 2.3.1 of [8] 

Definition 4.6. Two discretizations (L\,ip' 1 , a + ,a~), i = 1,2 have order r 
contact ifi> 2 (h,t,v) = i^ l {h,t,v) + 0{t r+l ) and L\{v) = L\{v) + 0(h r+1 ). 

Theorem 4.7. Let F^ be two discrete evolution maps of two discretizations 
(L l h ,tl) 1 , a + , a~), i = 1,2 of a regular Lagrangian system L: TQ — ► R. Then 
Ffo(v) — Ffr(v) + 0(h r+1 ) if ' {L l h ,ip l ,a + ,a~~) have order r contact. 



Proof. Assume the context and notations of the proof of Theorem 3.7 
mary, 

1 



t\h,v q )= { h 



Ll(v q ) 



>lh 



h^O, 
h = 0, 



are C k ~\ [/'CtxTQxTQis open, and 

C i = {(h,v,v)eU i :8»(v) = 8Jr(v)} 

are submanifolds of RxTQxTQ. Also, E and ( are defined, and (p 1 : C % 
by 



IxE 



0\h,v,v) 



h\c\d\+{i),d l -{v)) 



h,-(v + v, 



h^O, 
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7* : V' 1 — > J7 1 , where F l C R x £? is open, and J % (h, z q , —z q ) is the unique critical 
point in C % of L|<^ -1 (/i, z„, — z ? ). The maps F l are constructed as the graphs of 
7^307* where 7r 2 3 = (7r 2 , 7r 3 ) and 7r 2 , 7r 3 : IxTQxTQ — > TQ are the projections. 
Proposition 6 of j4 implies maps agree to the same order as their graphs, and 
one order higher if their residuals are symmetric. So it suffices to show that 
7T23 o 7 2 = 7T23 o 7 1 + 0{h r ) and that res r (7T23 o 7 , 7r23 o 7 1 ) is symmetric. 

To establish the contact order of 7r 2 3 ° 7 1 , the basic data of the corresponding 
critical point problems has to be compared. For that it is inconvenient that the 
manifolds C l depend on i. Let Q l be maps from R x TQ x TQ to itself which 
have the following properties: 

1. e 2 = e 1 + o(h r+1 ) ; 

2. 0\ i = 1, 2 are the identity on {0} x TQ x TQ; 

3. O', i = 1,2 have nonsingular derivatives on {0} x TQ x TQ; 

4. tq o 7r 2 o Q z (h,v,v) — d l ^~{v) and tq o 7r 3 o @ % (h,v,v) — 9^(v). 

For example, we can use a metric on Q the parallel transport P q2 . qi along 
geodesies between nearby points of Q x Q to define 

The purpose of the maps <d l is to normalize the submanifolds C % . In particular, 
by transporting the vectors defining the curve segments to common connection 
points, each 6* maps C % diffcomorphically to an open submanifold of R x TQ © 
TQ where 

TQ(STQ={ (v, v)£TQxTQ: tq{v) = tq(v) } 

is the Whitney direct sum. 

Since O l (0, v, v) — (0, v, v) and Q l is a local diffeomorphism on {0} x TQ x 
TQ, Theorem 1 of [4 implies that O l can be assumed to be a diffeomorphism 
between neighborhoods U\ U l D {0} x TQ x TQ. Set 

s i = s i o(e i )-\ <p l = $* o (e*)-\ f = e i of, 

where S(h,v,v) — L(h,v) + L(h,v). j(h, z q , —z q ) is the unique critical point 
of 5*|(R x (TQ © TQ)) in XJ l subject to the constraint (p = (h,z q ,-z q ). By 
the c onsis tency of L l h , the cons isten cy of ip 1 , the blow-up constructions of The- 
orem 3.7 and by Proposition 4.5 S2 = S\ + 0(h r ) and (p = <p + 0(h r ), so 
6 2 o 7^= 6 1 o 7 1 + 0(h r ). Also, since 6 2 = 6 1 + 0(h r+1 ), res r ( A f 2 ,j 1 ) = 
res r (7 2 ,7 1 ). So it is sufficient to show that res r (7T23 o 7 2 , ^23 o 7 1 ) is symmetric 
i.e. that 

r7r 2 res r (7 2 ,7 1 ) = Tvrares^^ 2 ,^ 1 ). 
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From Remark 3.8 setting h = 0, O l o 7*(0, z q ,— z q ) is the solution of the 



variational problem of finding the critical points of L(v) + L(v) with the con- 
straints v,v € T q Q and ^(v + v) — z q . This variational problem admits the 
Z 2 symmetry (v, v) i— > (v, v) and the solution is <d l o 7 J (0, z 9 , — z q ) = (0, z 9 , z q ) 
(and therefore all solutions occur on the fixed point set of the Z 2 action) . So it 
suffices to show that res r (S 2 ,S 1 ) and res r ((p 2 ,(p 1 ) are symmetric i.e. 



res r (S 2 , 5^(0, v, v) = res r (S* 2 , 5^(0, v, v), 
res r (^ 2 , ^(O, v, v) = res r (<^ 2 , ^)(0, t>, w). 



(4.3) 



The first of (4.3) is immediate: res r (S 2 ,S 1 ) is symmetric since S 1 and S 2 are, 
and 

res r (S 2 , S 1 ) = res'(5 2 o 6 2 , S 1 o 6 1 ) 

= res r (S 2 , S 1 ) oO^ TS 1 rcs r (6 2 , 6 1 ) 

= res r (S 2 ,S 1 )oe 1 . 



The proof of the second of (4.3 1 begins with the observation that 



(1,tq,tq) oQ^h^v) = (h,d^{v),di'(v)), 
so that, after defining the involution 

a(v,v) = (v,v), 
we have 

(h,di l + (t),di-(v)) = (i 7 T Q ,T Q )oO*o(i,<j)(h,v,v), 

and hence (by abuse of notation iT23{h, g/ + , q~) = (q + , <7 - )) 

7T23 ° <P* = T"C _1 ° ( T Q ° T2i T Q ° Ts) ° @* ° (lj Cr) o (6 l ) _1 . 



Using Proposition |4. 5 1 
vertres r (7r 2 o <^ 2 ,7r 2 3 o yj )(0,u, u) 

= vert res r+1 (C _1 o (r Q o 7r 2 , t q o 7r 3 ) o 6 2 o (1, cr) o (9 2 )" 1 , 

C 1 o (r Q o 7T2, r Q o 7T 3 ) o 6 1 o (1, CT ) o (6 1 )- 1 ) (0, w, S) 
= vert TC 1 T(T Q o 7 r 2 ,T S o^ 3 )(res r+1 (e 2 , 0^(0, 5,v) 

- T(l, a) res r+1 (0 2 , 6^(0, v, vj) . 



To compute the outer part of this, note that, for Sq + and Sq in the same fiber 
of TQ, 



veitTC- l (Sq + ,5q-) = ^(<Sc/ + - 6q~, -5q+ + Sq-) 
= 1 -(l-<j)(Sq + ,Sq-), 



(4.4) 
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and since 

a o T(t q o tt 2 , t q o 7r 3 ) (res r+1 (6 2 , 6 1 )(0, v, v) 

-T(l,a)res r+1 (6 2 ,e 1 )(0,w,w)) 
= T(t q o 7t 2 , t q o 7r 8 )T(l, a) (res r+1 (e 2 , G 1 )^, fi, «) 

- T(l, ff ) res r+1 (6 2 , G 1 )^, v, «)) 
= T(r Q o 7T 2 , r s o 7r 3 )(r(l, a) res r+1 (e 2 , 6 1 )(0, S, «) 

-res^+^G 2 ^ 1 )^, v,v)), 



it follows from (4.4 1 that 
vertres r (7T2 o <p 2 ,7T3 o (p )(0,v,v) 

= 2 T ( TQ 07r 2' T Q 07r 3)( 

res r+1 (6 2 , G 1 )^, v, v) + res r+1 (6 2 , G 1 )^, v, v) 

- T(l, a) res r+1 (e 2 , e 1 )^, v, v) - T(l, a) res r+1 (e 2 , 0^(0, 5, »)) , 

which is symmetric i.e. 

vertres r (7T2 o <p 2 , 7r3 o ip ){0,v,v) = vert res r (tt2 o v? 2 ,7T3 o ^ 1 )(0,w, v). 

This implies that res r (7T2 o (p 1 ,itz o (^ 1 )(0,w,w) is symmetric, since by Propo- 
sition |4.5| that is vertical anyway, so equality of the vertical parts is sufficient 
for equality. Thus, res r ((p 2 , (p 1 ) = (0,res r (7T2 o (p 2 ,iT3 o (p 1 )^ is symmetric, as 
required. D 

Theorem 4.8. Suppose that Fh be an evolution map of an order r discretiza- 
tion (X/j, ip, a + ,a~) of a regular Lagrangian system L : TQ — ► R. Then Fh(v) = 

F^(v) + 0(h^). 



Proof. Combine Theorems 4.3 and 4.7 □ 



5 Qx Q 

In this Section we provide the relations between the discrete mechanics of Defi- 
nition [33J with discrete phase space TQ, and the standard discrete mechanics, 
with discrete phase space Q x Q. By standard discrete mechanics we mean: 

1. The discrete Lagrangian is of the form L^ : Qx Q — > R. 

2. (h,(q + ,q~),(q + ,q~)) is critical if ((g + , g~), (g + , g - )) is a critical point of 
the variational principle 

f ^,? XQ ((9 + , «"), (t,r)) = L^(q+,q-) + L%**{q+, q'), 
I q~ and q + constant, and q + = q~ . 
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3. The discrete evolution is defined to advance from (q + , q ) to (q + , q ) such 
that (h,(q + ,q~),(q + ,q~)) is critical. 

4. Given a continuous Lagrangian system L : TQ — > K, a discrete Lagrangian 
is order r if 

L(v) = L^ q (t q o F**(v),t Q (v)) + 0(h r+1 ), 

in which, of course, any order r accurate approximation F h E + 0(h r+1 ) 
may be substituted for the exact flow F h E . 

This is the variational principle used in the example context of Section [2] and 
also in [8]. 

The equivalence of the two formalisms is based on Lemma |5.1[ which is 
obtained in the proof of Proposition 2.9 of [I]. 

Lemma 5.1. Let M. be a manifold and let (ip,a + ,a~) be a discretization of 
the tangent bundle of M. Then \&(h, v) = (h, d£(v), d^(v)) is a diffeomorphism 
between open neighborhoods U \ ({0} x TM) and W \ ({0} x (M x M)). 

Given L^ xS , choose any order r discretization (tp,a + ,a~) of Q, define 

L h (v)^Lf* Q (d+(v),d-(v)), (5.1) 

and consider the discrete Lagrangian system (Lh,ip,Ci + , a~). This is order r by 
Definition |4.6[ and ^ intertwines the objectives and constraints of the discrete 
variational principles on TQ and Q x Q. Therefore ^ is a bijective correspon- 
dence between their critical points. Conversely, if (Lh,tpi a+ : a ~) is a discrete 
Lagrangian system, then define L^ x such that ( |5.1| holds, and the same cor- 
respondence is obtained. Thus, the version of discrete mechanics on TQ and 
the standard version on Q x Q arc entirely equivalent, in an order-preserving 
way. 
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